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We consider scalar-tensor theories and reconstruct their potential (7('I>) and coupling -F($) by 
demanding a background ACDM cosmology. In particular we impose a background cosmic history 
H{z) provided by the usual flat ACDM parameterization through the radiation {weff ~ 1/3), 
matter (weff = 0) and deSitter (we// = —1) eras. The cosmological dynamical system which 
is constrained to obey the ACDM cosmic history presents five critical points in each era, one of 
which corresponding to the standard General Relativity (GR). In the cases that differ from GR, 
the reconstructed coupling and potential are of the form F(<I>) ~ <I>'^ and (7($) ~ _F(<I?)'" where m 
is a constant. This class of scalar tensor theories is also theoretically motivated by a completely 
independent approach: imposing maximal Noether symmetry on the scalar-tensor Lagrangian. This 
approach provides independently: i) the form of the coupling and the potential as F{^) ~ $^ and 
[/($) ~ F($)™', ii) a conserved charge related to the potential and the coupling and Hi) allows the 
derivation of exact solutions by first integrals of motion. 

PACS numbers: 98.80.Es,98.65.Dx,98.62.Sb 
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I. INTRODUCTION 

There is accumulating observational evidence based 
mainly on Type la supernovae standard candles [l| and 
also on standard rulers 0, Q that the universe has en- 
tered a phase of accelerating expansion at a recent cos- 
mological timescale. Such an expansion implies the exis- 
tence of a repulsive pressure on cosmological scales which 
counterbalances the gravitational attraction of matter on 
these scales giving rise to an overall accelerating behav- 
ior. There have been several theoretical approaches (see 
0, [E] for a review) to understand this phenomenon. The 
simplest of such approaches assumes the existence of a 
positive cosmological constant which is small enough to 
have started dominating the universe at recent times. 
The predicted cosmic expansion history, in this case, (as- 
suming flatness) is 



H{z) 



where fior = ^737 



(1.1) 

10 is the energy density of radi- 
ation today normalized over the critical density for flat- 
ness Per it- Also Oom = — 0.3 is the normalized 
present matter density and = 1 — fiom ~ is the 
normalized energy density due to the cosmological con- 
stant. This model provides an excellent fit to the cos- 
mological observational data and has the additional 
bonus of simplicity due to a single free parameter. De- 
spite its simplicity and good fit to the data, such a model 
fails to explain why the cosmological constant, so unnatu- 
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rally small dominates the universe at recent cosmological 
times. This problem is known as the coincidence prob- 
lem. Furthermore, there is no urgent theoretical reason 
implying IIa — 0.7 and f2om — 0.3 at present time so also 
a fine tuning problem has to be taken into account. 

In the effort to address these problems, several mod- 
els, which essentially can be grouped in two classes, have 
been proposed: The first class assumes that General Rel- 
ativity (GR) is valid at cosmological scales and attributes 
the accelerating expansion to a dark energy component 
which has repulsive gravitational properties due to its 
negative pressure. The role of dark energy is usually 
played by a scalar field minimally coupled to gravity 
called quintessence Alternatively, the role of dark 
energy can be played by various perfect fluids (eg Chap- 
lygin gas 0) topological defects [1], holographic dark 
energy etc. 

The second class of models attributes the accelerating 
expansion to modifications and extensions of GR which 
converts gravity to a repulsive interaction at late times 
and on cosmological scales. Examples of this class of 
models include scalar-tensor theories [l3, [IH, fiM) 
tended gravity theories [l3], braneworld models [ll] etc. 
An advantage of models in this class is that they natu- 
rally allow [3, [l5| a super-accelerating expansion of the 
universe where the effective dark energy equation of state 
w — ^ crosses the phantom divide line w = —1. Such 
a crossing is consistent with current cosmological data 

A representative model of the second class is provided 
by scalar-tensor theories of gravity. In these theories, 
the Newton constant is obtained by dynamical properties 
expressed through the coupling F(<^). Dynamics is given 
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by the Lagrangian density [H, [11, [U 

, (1.2) 

where C^[ipai'i dfj.!^] represents matter fields approximated 
by a pressureless perfect fluid in the dust dominated 
regime. We have set SttG = 1 and e — ±1 for standard 
scalar and phantom fields respectively, i.e. we have nega- 
tive kinetic energy for the scalar degree of freedom. Note, 
however that negative energy needs e = — 1, but there 
also exist positive energy configurations with e — 
This happens due to the fact that the kinetic term of the 
actual spin-0 degree of freedom does not come only from 
the obvious (9^<1>)^, but also from the cross term F($)i?, 
since R involves second derivatives (see [2l| for details). 
As discussed below, for e = and <I> ^ i? the Lagrangian 
(|1.2p can also describe f{R) generalizations of GR. 

In the present study, we reconstruct the potential C/(<&) 
and the coupling F{^). Instead of specifying various 
forms of U{^) and F{^) and finding the corresponding 
cosmological dynamics, we specify the cosmological dy- 
namics to that of the ACDM cosmology and search for 
possible corresponding forms of J7($) and F{^). The 
original method for the reconstruction of scalar-tensor 
theories from a given cosmic history H(z) was introduced 
in Ref. [l^ and applied to specific cases of late cosmic 
history in Ref. [ij, [2l| . Our reconstruction approach is 
different in two aspects: 

• We use a dynamical system formalism and find the 
critical points that determine the generic evolution 
of the system. 

• We start the reconstruction from the radiation era 
rather than focusing only at late times through the 
acceleration epoch. 

In particular, we consider the general dynamical system 
for scalar-tensor theories and study the dynamics oiU{^) 
and F($) using as input a ACDM cosmic expansion his- 
tory. Our study is performed both analytically (using 
the critical points and their stability) and numerically 
by explicitly solving the dynamical system. 

Since our reconstruction assumes a fixed cosmic his- 
tory background the stability of the critical points we 
find should be interpreted with care. The fixing of the 
cosmic history has eliminated perturbations of H{z) but 
has allowed for perturbations in the forms of the coupling 
F($) and the potential t/($). In the usual stability ap- 
proach the coupling F($) and potential C/($) are fixed 
and perturbations are allowed in the cosmic history H{z) 
to determine the stability of the phase space trajectories. 
Since the later stability approach is more physical (but 
docs not lead to a reconstruction) we focus on the actual 
values of the critical points and interpret their stability 
only as a test of the corresponding numerical evolution of 
the same system. Thus, all the critical points we find are 
assumed equally important cosmologically, independent 
of their stability. 



The structure of the paper is the following: In the next 
section we derive the dynamical system for the cosmolog- 
ical dynamics of scalar-tensor theories. Using as input a 
particular cosmic history H{z) (eg ACDM), we show how 
this system can be transformed so that its solution pro- 
vides the dynamics and the functional form of U{^) and 
F($). We also study the dynamics of this transformed 
system analytically by deriving its critical points during 
the three eras of the cosmic background history (radia- 
tion, matter and deSitter). We find that the cosmological 
dynamical system constrained to obey the ACDM cosmic 
history has five critical points in each era, one of which 
corresponds to GR. In Sect. Ill, we use the solution of 
the above system to reconstruct the cosmological evolu- 
tion and functional form of the coupling and potential, 
which are of the form F($) $^ and U{^) ^ F($)'" 
where m is an arbitrary constant. We show that such 
forms are also motivated by a completely independent 
approach ie by imposing Noether symmetry on the scalar- 
tensor Lagrangian We also demonstrate the agree- 
ment between the analytical and numerical results of our 
reconstruction scheme. Finally, in Sect. IV, we conclude, 
summarize and refer to future prospects of this work. 

II. DYNAMICS OF SCALAR-TENSOR 
COSMOLOGIES 

Let us consider the action (|1.2p describing the dynam- 
ics of Scalar Tensor theories in the Jordan frame. In the 
context of flat Friedman-Robertson- Walker (FRW) uni- 
verses, the metric is homogeneous and isotropic i.e. 

ds^ = -dt^ -f a2(t)dx2 (2.1) 

Variation of the action (|1.2p with respect to the metric 
leads to the following dynamical equations which are the 
generalized Friedman equations [20, 2^, [23j 

3FH^ = pn, + Pr + ^e'i>^ -3HF + U (2.2) 

-2FH = p„, + ^pr + €<i>^+F-HF (2.3) 

and variation with respect to $ gives the Klein-Gordon 
equation: 

e($ + 3i/$) = 3F($),$(ij + 2H'^) - [/($),$ (2.4) 

where we have assumed the presence of perfect fluids pm , 
Pr representing the matter and radiation energy densities 
which are conserved according to 

An + SHpn, = , (2.5) 

Pr + 4:Hpr = . (2.6) 

The background equations in eqs. (|2.2p and (|2.3p can 
be rewritten in a more convenient form, which is easier 
to confront with observations (see for example Ref. 24]): 
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BFoF^ = pDE + Pm + Pr, (2.7) 

4 

- 2FaH = + pde + Pm + ^Pr (2.8) 
where we have set 

PDE = \ef!?'^ -'iH'^{F-Fo)-ZHF + U (2.9) 



where we can associate a;4 with and Xi + X2 + £2:3 = 
f^DB with curvature dark energy (dark gravity). Defining 



also Vlr. 



3FH 



2 we can write equation (|2.14p as 



n„i = 1 — xi — X2 — ix\ — X4 (2.20) 



Let us now use (|2.15p to express (12.31) as 

F,'2 



ll 



2 pr F" H'F' F' 



2FiJ2 3FiJ2 -'2F 2F 2FH ' 2F 

(2.21) 



or 



PDE = -^■' + F + 2HF-U-{2H+iH'%F,^F) (2.10) ^ 3_2xi-3x2 + X4 + 3£x^ + x? + 2:^~xi- (2.22) 



H' 



and the subscript "0" denotes present day values. The 
function pjj e defined in this way can be shown to satisfy 
the usual energy conservation equation: 

PDE + iH{pDE + pde) = (2.11) 
where the DE equation of state is defined as 



WDE 



PDE 



-1 



F - HF - 2H{Fo - F) 



PDE ^e<^>^ -SH^{F - Fo) -3HF + U 

(2.12) 

By using eqs. (|2.7p and (|2.8p we can express the equa- 
tion of state wde 3.S 



wde = 



3E{z) - {\ + z){dE{z)/dz) 
3i;(z) -3r!l?(l + z)3 



(2.13) 



where we have set E{z) = H^{z)/Hq. This relation is 
exactly the same as in standard Einstein gravity Q and 
therefore we can constrain wde from Type la supernovae 
observations in the same way. 

In order to study the cosmological dynamics implied 
by equations (|2.2p . (|2.3p and (|2.4p we express them as a 
dynamical system of first order differential equations. To 
achieve this, let us first write ()2.2p in dimcnsionless form 
as 



Pr 



U 



3FH^ 3FH^ 6F 3FH^ 



El. 



(2.14) 



where 



(2.15) 



, d _ d Id 
^ dlna ^ dN ^ H'dt 

Let us now define the dimcnsionless variables xi,...,Xi 
as 



Xl 
X2 

X4 



EL 

'~F 

u 



3FH^ ' 
$'2 

Pi- 
3FiJ2 



(2.16) 
(2.17) 

(2.18) 
(2.19) 



H H 

Differentiating X4 of (|2.19p with respect to N , we have 

P'r Pr F' 2pr H' 



3FH^ 3FH^ F 3FH^ H 



x'a = —4x4 + a;4a;i — 2x4 



ll 



(2.23) 



(2.24) 



where we have used (|2.6p . Similarly, differentiating (|2.17p 
with respect to N, we find 



X2 



where 



Xl (1 — m) 



F^^/F 



H' 
) 

' H 



(2.25) 



(2.26) 



and ,$ implies derivative with respect to $. Finally dif- 
ferentiating (|2.18p with respect to N and using (|2.4p . one 
finds 



.H' 



H' 



2exn Xl — 



H 



H 



e.{x\y — ex\xi — &f^x\ — 2xi -f mx2Xi 

(2.27) 

The dynamical system ^(TM . (I^TM)) . ^T^, (P?^ de- 
scribes the cosmological dynamics of Scalar- Tensor theo- 
ries. In two special limits, it can be transformed into 
the dynamical systems obtained in /(i?) theories and 
Quintessence respectively [l^. However, we have to say 
that there is one more degree freedom in scalar-tensor 
theories, compared to /(i?) theories, since here we have 
two arbitrary functions, i.e. F{^) and U(^). In general, 
we should add a further parameter n which could be re- 
lated to Fq, and -F',*^, beside the above m. In fact, if 
we do not try to reconstruct the function F{<^)^ such a 
function can be fixed a priori and the corresponding pa- 
rameter n would be, for example, F^^ / F ^ or some similar 
condition. In such a case H{N) would not be fixed as in 
our reconstruction approach but would have to be deter- 
mined by the autonomous system. This is the approach 
followed in [26| where F^/F is not present as a variable 
in the autonomous system since F{^) is fixed a priori. 



4 



The parameters in [2q are the power law of the fixed po- 
tential (there called n) and ^, the couphng F(<i>) — ^<i>^. 
In that case, H{N) and its perturbations are allowed to 
vary. However, in a reconstruction approach where H{N) 
is fixed a priori, while F_^/F (and its perturbations) are 
allowed to vary, -F($) can be reconstructed from the au- 
tonomous system. Thus, in our reconstruction approach 
where H(N) is fixed, n is no more a parameter and it 
is allowed to evolve. In this sense, the present approach 
is more general than that presented in [2^. Finally, as 
we shall see below, the reconstructed function F{^) and 
[/(<&) can be independently obtained asking for a sort of 
" first principle" , the existence of a Noether Symmetry. 
Since the results coincide (the form of F and U are the 
same in both the approaches), we are confident that the 
presented method is consistent. 

Considering again the /(i?)-theories, by setting U = 



[27| and using the transformations 



Xi - 


Xi 


(2.28) 


X2 - 




(2.29) 


X4 - 


-> X4 


(2.30) 


H' 




(2.31) 


H 


-y X3~2 


$ - 


R 


(2.32) 


e - 


^ 


(2.33) 



one recovers the dynamical system of f{R) theories (eqs 
(2.15), (2.17) and (2.21) in Ref. ^ ), where the tilded 
{') quantities are the ones defined for f{R) theories: 



and A = -U^^/U. 

It is worth noting that the results of our analysis do 
not rely on the use of any particular form of H{z). They 
only require that the universe goes through the radiation 
era (high redshifts), matter era (intermediate redshifts) 
and acceleration era (low redshifts). The corresponding 
total effective equation of state 



eff 



-1 



2H'{N) 
3 H{N) 



(2.43) 



IS 



"eff 



— Radiation Era 
o 



ygjy — Matter Era 
Ueff = —1 deSitter Era 



(2.44) 



For the sake of definiteness however, we will assume a 
specific form for H{z) corresponding to a ACDM cos- 
mology (|1.1[) which, in terms of N , takes the form 

H{Nf = Hi [riome-^^ + riore"^^ + ^Ja] (2.45) 

where N = Ina = — ln(l + z) and = 1 — fiom — ^^Or- 
Also, we can use (j2.45l) to find ^h{n) ' ^ quantity needed 
in the dynamical system, as 

H'{N) _ -31}o„e-3^-417o^e'-^^ 

H{N) " 2 (1 - 17o™ - ^or + + nare-^N) 

(2.46) 

The crucial generic properties of are their values 

at the radiation, matter and deSitter eras: 



T — is — 3x2 + xl+ X4 



(2.34) 



-i^j - 2^3(^3 - 2) - X2(2i3 - ii - 4) (2.35) 



-2X3 Xi + Xi X4 



(2.36) 



On the other hand, the following set of transformations 
gives the autonomous system for Quintessence (see eqs 
(175) and (176) in Ref. @): 



Xi 
X2 



X4, 



X 





with dynamical equations 
x' = -3x+^eXy'^ - 



-x(ex 



y' = --^Axy -I- ^y{ex^ + l-y^) 



(2.37) 
(2.38) 
(2.39) 
(2.40) 



(2.41) 



(2.42) 



H'{N) 
H{N) 
H'jN) 
H{N) 
H'jN) 



= -2 iV < Nr. 



(2.47) 



-- N„n <N<N^A (2.48) 



= N>N,nA 



(2.49) 



where Nr 



-In 



and N„ 



,™ _ -...-f^ -..^ ^v,„A "il^n^ are the 
N values for the radiation-matter and matter-deSitter 
transitions. For flom = 0.3, fior = 10~^ we have Nrm — 
—8, A^mA — —0.3. The transition between these eras 
is model dependent but rapid and it will not play an 
important role in our analysis. 

It is straightforward to study the dynamics of the sys- 
tem (I2:22l) . ((2:24)) . ((2:251) . by finding the criti- 
cal points and their stability in each one of the three 
eras. Notice that even though this dynamical system is 
not autonomous at all times, it can be approximated as 
such during the radiation, matter and deSitter eras when 



H'jN) 
H{N) 



is approximately constant. The critical points and 
their eigenvalues are shown in Table I. An interesting fea- 
ture to observe in Table I is that in each era there are 
five critical points but only one of them is a stable "at- 
tractor" for a given value of m. Also, remembering that 
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TABLE I: The critical points of the system (|2.22|) . (|2.24l) . (I2.25|l . (|2.27|l and their eigenvalues in each one of the three eras (rad. 
era iV < -In^^, matter era -In^^ < iV < -ilnTp^, deSitter era N > -iln^r^). 



H/ra 






22 


a.2 
2^3 


2:4 







Eigenvalues 




Ri 


2 





-1 








1 


(2,3,l,6-2m) 




R2 


1 














1 


(1,2,-1,5-m) 


Radiation 


Rs 


-1 











2 


-1 


(-1,-2,-3,3-hm) 


1 

■Weff = 3 


Ri 


4 


15~8m-|-m^ 


2(771 — 5) 777. 








1 


V m — 1 ' m — 1 ' / 


-1 + m 


3{m-l)^ 


3(m-l)^ 




Us 


U 


U 


U 


1 
i 


A 

u 


n 

u 


11,-1,-2,4] 




Ml 


2 





-1 








1 


(1,2,1/2,5 -2m) 




Ma 


3/2 





-1/2 








1 


(l/2,3/2,-l/2,-3/2(m-3)) 


Matter 


Ma 




















(-1,-3/2,-2,3) 


U)e/7 = 


Mi 


3 


15 — llm + 2m^ 


1_ 9^71 -|_ 2777,^ 








1 


/ 4- — 777 M h\ 

(^,;;i^,see ') 


m — 1 


4{m-l)^ 


4(m-l)^ 




Ms 


1 





-1/4 


1/4 





3/4 


(1,-1/2,-1,4 -m) 




Ai 


2 





-1 








1 


(-2,-1,-1,2 -2m) 




Aa 


3 





-2 








1 


(-1,0,1,3 -3m) 


deSitter 


A3^ 


3 





-2 








1 


(-1,0,1,3 -3m) 


Weff = -1 


A4 





1 











1 1 


f A V24m+l + 5 v'24m+l-5\ 
;^ ' 2 ' 2 ^ 




As 


4 





-4 


1 








(1,1,2,4 -4m) 



-63m2+llgm+l-ll — 3777. + Y 8m^ -63m2 + 118m + l + ll 

" the other two eigenvalues are: >^ tt, r; , n n 

^ 2(m — 1) ' 2m — 2 

' the other two eigenvalues are: - 7^+^48^3-2^^2+358^+1-19 _ -7m+V48m3 -2^^2+358^+1+19 
Notice that A2 and A3 are degenerate 



the positivity of the energy of the (helicity zero) scalar 
partner of the graviton, i.e. the positivity of the kinetic 
energy of the scalar field (see Ref. [lH) is expressed by: 

then we have in dimensionless variables 



-f+xl>Q (2.51) 

All points of Table I satisfy (|2.5ip except i?4 and M^. 
In order for i?4 to satisfy (|2.5ip . we must have m < 2 or 
m > 3 while for M4 it is m < 2 or m > 5/2. Therefore 

the allowed region for to, so that we have critical points 

2 

that are physical + x\ > 0), is 

TO < 2 or TO > 3 (2.52) 

Now, regarding the "attractor" behavior of the system 
in each era separately, we see that in: 



• Radiation Era: i?3 is an "attractor" for to < —3 
while i?4 is an "attractor" for —3 < m < 1; 

• Matter Era: Af4 is an "attractor" for to < 1; 

• deSitter Era: Ai is an "attractor" for to > 1 and 
A4 is an "attractor" for to < 1. 

As discussed in the introduction, the significance of the 
"attracting" nature of the critical points is limited be- 
cause the reconstruction method we are using has forced 
us to eliminate physical perturbations of H{N) and has 
introduced unphysical perturbations of the functional 
forms of the potential C/($) and coupling F{^). How- 
ever, we can use the "attractor" critical points as a pre- 
diction for the numerical evolution of our reconstruction 
dynamical system. 

A particularly interesting feature of the critical points 
of Table I is that in all the cases that differ from GR, 
the expansion rate in each era is induced by dark energy 
{^DE — 1) which implies that the scalar field $ could 
also play the role of dark matter if it is found to have 
the proper perturbation properties at early times. We 
postpone the analysis of such perturbations for a future 
study. 
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0.2 


-0.2 

g 

i -0.4 
-0.6 
-0.8 
-1 



-10 

N 



FIG. 1: The effective equation of state ■Weff(N) imposed on 
the dynamical system (obtained from (|2.43p using (|2.45p '). 



$(7V) = -2V6 ^ i^o/'e-^i^/2 + C 

Xi 
X3 pl/2 



= 2V6f Fl'^ (l - e-*^^/2 j + $0 (3.2) 



where 



Xi 



1 xl 



(3.3) 



and $0 = ^'(^ = 0). Equations ([3TT|l and ([3?2|l allow 
us to eliminate N in favour of $ 



(3.4) 



To confirm the dynamical evolution implied by the "at- 
tractors" of Table I, we have performed a numerical anal- 
ysis of the dynamical system (|Z!22)) . ((2?24)) . ((2?25)) and 

(P:?f|) using the ansatz (PIIS]) for with f^om = 0.3 



leads to the 



and r^or = 10 ^. This ansatz, for ^jjj^, 
Weff{N) shown in Fig. 1. We have set up the system 
initially, on i?4, with m = —0.5. As seen in Fig. 2 
and Fig. 3 the system follows the evolution of R4 from 
the Radiation Era, to the Matter Era (M4) and finally 
to the deSitter Era (A4). We have checked that if we 
choose initial conditions not exactly coinciding with any 
of the other critical points then the system is captured by 
the i?4 "attractor" and follows the trajectory mentioned 
above, i.e. (init) R4 M4 — > A4. 

Finally, if the initial conditions of the system are set 
to be exactly on the top of other critical points then the 
system will end up on the Ai critical point. 

The choice of a constant m is justified by the fact that 
Xl and a;| are constants in each era (see Fig. 2). Also, the 
potential U{^) and coupling F{^) that are mostly used 
in the literature are power-laws or exponentials, which 
also give a constant m (see equation (|2.26p ). This feature 
will be obtained also by Noether symmetries, as we will 
discuss below. 



III. RECONSTRUCTION OF F($) AND i7($) 

A. Analytical Results 

Our goal is now the reconstruction of the form of the 
potential U{^) and coupling F($) corresponding to each 
one of the critical points of the system shown in Table I. 
Let us consider a critical point of the form (ii, X2, ^3, 2:4). 



Using (|2T6l) . we find 



F = Fne^^i^ 



(3.1) 



where Fq = F{N = 0) is the present value of F. Using 
now eq. (|2.18p we find 



where = ^fi- The quadratic form of F($) can be 
achieved in a completely different approach imposing 
Noether symmetry in the scalar-tensor Lagrangian. 
From eq. (|2.17p . we have 



U{N) = X2 ■ 3F{N)H{Nf 



(3.5) 



Using now the input form of H{N) (Eq. (|2.45|) ). we find 
the dominant term of H{N) in each era, that is 

( l^or-e""^, Rad. Era 

H{NflHl = <^ l^ome"^^, Mat. Era 

[ 1 — f^or — i^Om, dS Era 
Thus using eq. ([X^ we have ff($): 

' ^[i{^-CyY'"\ Rad. Era 



Mat. Era 

1 — r^or — f^m, dS Era 

and then we can express (|3.5[) in terms of $ to get the 
relevant form 



C/($) = A($ - C) 



2+a 



(3.6) 



where, 
A= < 



f 3a;2^C'+^/^% 



Rad. Era 
Mat. Era 



, 3x2^(1 — f^or — f^om), dS Era 

{8/xi, Rad. Era 
6/xi, Mat. Era 
0, dS Era 
By using eqs. (|2.26p . p.4p and (|3.6p . it is straightfor- 
ward to find 



2 + a = 2m 
so that Eq. p.Gp can be written as 

[/($) = A($ - C)2™ 



(3.7) 



(3.8) 



Notice that even though the eigenvalue xi changes in the 
sequence i?4 — > M4, the exponent 2+a remains constant. 
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FIG. 3: The phase space trajectories on the xi — X2 plane (FigSa) and X3 — X4 plane (FigSb). 
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Furthermore, since m is constant (see Eq. p.7p ) 
Ea. (|2.26p allows to write U in terms of F, i.e. 



(3.9) 



The above analysis is only valid when the parameters 
xi, X2 and x§ are not equal to zero, as for example is 
the case for i?4 or when we have perturbed the initial 
conditions around a critical point. Otherwise, we have 
the following cases: 

• X2 = S§ = 0, then $ = $0 ^ const, F = F{N) 
given from eq. (|3.ip and U = 0, as is the case for 
i?2 and i?3. 

• Xi = = and we are in the deSitter era, then 
$ = $0 = const, F = Fq and U = Uq, as is the case 
for A4. 

9 xi = X2 = x'^ — 0, then $ = (f>o —const, F = Fq 
and U — 0, as is the case for i?5. 

Note that if we were considering other (non-critical) 
points, not only the reconstruction would involve a time 
dependence, but we would not be able to reconstruct si- 
multaneously F{^) and U{^) from the single function 
p.ip . Two arbitrary functions obviously need two ob- 
served functions to be reconstructed in the general case 
(see [3, [13, mi ) ■ In our case, however, the fact that the 
reconstruction occurs on the critical points means that m 
is fixed to a constant and this closes the system of equa- 
tions (11211), (ESI), (H^Sl), allowing us to proceed 
with the reconstruction numerically and analytically. 

Also, it is easy to see that the reconstructed theo- 
ries are merely Brans-Dicke theories with an additional 
potential U. If we define F = f](j)BD — then 

$ — J j4>BD and the Lagrangian (|1.2p becomes 



C 



•)BD 



1 t^BD 

2 



dbd 



>BDd„(f)BD ~ U {<j)BD) + Cm 

(3.10) 

where lubd = = constant. 

A noteworthy feature of the above reconstruction 
scheme is that Eq. p.9p along with Eq. p.4p are exactly 
the conditions for the existence of a Noether symmetry 
in Scalar- Tensor theories as we discuss in what follows. 



B. Noether Symmetries in Scalar- Tensor Gravity 

1. Generalities on the method 

Solutions for the Lagrangian (|1.2p can be searched by 
the so called Noether Symmetry Approach . This ap- 
proach allows, in principle, to find out cyclic variables 
related to conserved quantities and then to reduce dy- 
namics. Besides, the existence of symmetries fixes the 
forms of the coupling F{^), of the potential [/(<&) and 
gives the relation between them. 



Let us give a quick summary of the approach for finite 
dimensional dynamical systems before the application to 
our specific problem. 

Let C{q^,q^) be a canonical, non-degenerate point-like 
Lagrangian in the configuration coordinates g* (the "po- 
sitions"), where 



dC 

-—- = ; detHij = det 



^0. (3.11) 



Hij is the Hessian matrix related to C. The dot indi- 
cates derivatives with respect to the affine parameter A 
which, in general, corresponds to the time t. We are 
going to consider only transformations which are point- 
transformations. Any invertible and smooth transforma- 
tion of the "positions" Q' = Q' (q) induces a transforma- 
tion on the "velocities" such that 



Q'(q) 



(3.12) 



The matrix J' — \ \dQ'' / dq^ \ \ is the Jacobian of the trans- 
formation on the positions, and it is assumed to be 
nonzero. The Jacobian of the "induced" transforma- 
tion is easily derived and it has to he ^ Q ^ 0. 
Usually, this condition is not satisfied in the whole space 
but only in the neighbor of a point. It is a local trans- 
formation. A point transformation — Q*(q) can de- 
pend on one (or more than one) parameter. In general, 
an infinitesimal point transformation is represented by 
a generic vector field acting on the space The 
transformation induced by p.l2p is then represented by 

X is called the "complete lift" of X \2^. A function 
/(q, q) is invariant under the transformation X if 

ix/.o'wf +(^a'(q))f ^0, (3^4) 

where Lx/ is the Lie derivative of /. In particular, if 

Lx-C^O, (3.15) 

X is said to be a symmetry for the dynamics derived from 
the Lagrangian C. To see how Noether's theorem and 
cyclic variables are related, let us consider a Lagrangian 
C and the related Euler-Lagrange equations 



d dC _ dL_ 

dX dqi dqi 



0. 



(3.16) 



Let us consider also the vector field (|3.13p . By contract- 
ing (|3.16p with the a"s, one obtains 



d dC _ dC_ 

dX dqi dq^ 



= 0. 



(3.17) 
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Being 



d dC _ d_( jdC \ _ fda^\ dC_ 



dX dq3 
from p.l7p . we have 



(3.18) 



dX \ dq' 



(3.19) 



As a consequence, the Noether Theorem enunciates: 
If LyiC — 0, the function 



:dC_ 



(3.20) 



is a constant of motion. 

It is worth noting that Eq. (|3.20[) can be expressed in- 
dependently of coordinates as a contraction of X by a 
Cartan one-form 



3 - J » 



(3.21) 



Thus Eq. p.20p can be written as 

i^Oc - So . (3.22) 
where ix is defined through the relation 

ixrfg' = a' (3.23) 
By a point-transformation, the vector field X becomes 
d 



X - (zxrfQ") 



(3.24) 



X is still the lift of a vector field defined on the "space 
of positions" . If X is a symmetry and we choose a point 
transformation such that 



ixdg^ = 1 ; ixd<9* = i^l 



we get 



d 



dC 



(3.25) 



(3.26) 



Thus is a c ycli c coordinate and the dynamics can 
be reduced [l^, HH. Clearly the change of coordinates 
defined by (I3.25P is not unique. Usually a clever choice is 
very important. It is possible that more than one vector 
field X is found. In this case, more than one symmetry 
exists. 



2. The case of Scalar-Tensor Gravity 

The above method can be used to seek for solutions in 
the dynamics given by Lagrangian (|1.2p . In particular, 



for flat FRW metric, the field Lagrangian (|1.2p reduces 
to the point-like Lagrangian 



C = -3aa^F-3F$$a^d+a^ - U{^) ]-Da 



1 



-3(7-1) 



(3.27) 

where, for the sake of simplicity, we are considering only 
the scalar field case (the generalization to the phantom 
field case is obvious). The constant D is related to 
the perfect-fiuid matter density, being — D{ao/a)^'^ , 
where 1 < 7 < 2 defines the Zel'dovich range for the 
equation of state of standard matter. The above dy- 
namical system (|2.2[) , (|2.3p , (|2.4p is immediately deduced 
considering the energy condition and the Euler-Lagrange 
equations for p.27p . In the case of standard dust matter, 
7 = 1, the last term in (|3.27p reduces to an additive con- 
stant. Being {a, $} the configuration space of the system, 
the problem is 2D and then the infinitesimal generator 
of the Noether symmetry is 

oa o9 da 
where a and [3 are functions depending on a and and 



da , 
da 



da ■ 
— $ 



dl3_. 
da 



9$ ■ 



(3.29) 



The condition for the existence of a Noether symmetry 
is Ly^C — Q. It, explicitly, gives an expression of second 
degree in d and $, whose coefficients are zero due to the 
fact they are considered to be linearly independent. Then 
this set of coefficients gives rise to the following system 
of partial differential equations (2^ , 

a + 2a— + a- — + a/3— = (3.30) 

da oa I' t 



2. + 4 + «1) + <,./3 + 2^^ 4f = 
oa d9 J d9 3 oa 

(3.31) 

3a-6F$|| + 2a|^ = (3.32) 



3a 

IT ~ 

Equation (|3.33p can be re-written in the form 



(3.33) 



IT 



where 



m ■ — — 

F 

3a F 



(3.34) 



(3.35) 



It is worth noting that Eq. p.34p is a relation between 
the potential and the coupling and it exactly coincides 
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with (|2.26p . Solving the above system means to find out 
the explicit form of the set of functions {a, (3, F, U}. For 
this purpose we consider the separation of variables, 



Aiia)A2i'P) 



and 



(3.36) 



and F. Using eqs. ([O0| and (|04| for B2 and B'^ on 
eq. p.45p we get 



Also, using (P^ii)) in ((g^ yields 



(3.46) 



/? = Bi(a)B2($) 
Then from eq. p.32p we get 



Bia 



A2_ 
AB' 



A',F' 
2B'^ 



C 



(3.37) 



(3.38) 



where C is a separation constant. The solution of (|3.38p 
is simple and we get 



Ax 



Bia 
'6C 



(3.39) 



and 



A2 - 2A'2F' 
4C 



(3.40) 



Using eqs. ((3:36)) and ((3:37)) on ((330)) we get 



a dBi 3{A2F + 2CB2F') 



Bi da 2{A2F + 3CB2F') 



(3.41) 



where s is a separation constant. Hence, from eqs. (|3.39p 
and (I3.4ip we get 



Bi = Ba" 



and 



Ai 



B 

6C' 



,s+l 



(3.42) 



(3.43) 



where B is a constant of integration. Also, from (|3.4ip 
we have 



EL 
~F 



2s + 3 A2 



QC{s + 1) B2 
Using eqs. ((SHj), ((g:i5)) and ((^3^ yields 



(3.44) 



2FA'2 + {A2{s + 3) + QCB'2) F' - 2B2C {s - 3F") = 

(3.45) 

Now we are left with three equations (|3.40p , p.44p and 
(I3.45P to solve for the three unknown functions A2, B2 



2A'2F' (3(5 + 1) {F'f + F{23 + 3)) + 
+A2 ({s + 3) {F'f - 2F(2s + 3)F") = (3.47) 



Now, we can use equations p.46p and ()3.47p to elimi- 
nate A2 and A2 in favor of F, 



F" 



3s(s + l)(s + 2)F"* 

(2s + 3)F2 
(s + 1) (8s2 + 16s 



i)F'^ s(2s + 3^ 



2(2s + 3)i^ 



■3.48) 



Equation (|3.48p is nonlinear and its complete solution 
is an elliptical integral of second kind which is not simple 
to handle. However, an exact solution can be found to 
be of the form 



(3.49) 



Using the ansatz (lOO)) in ((OS)) we get that ^ = 
or ^ — — i, with the latter corresponding 



24(s+l)(s+2) 



to the conformal coupling [Sj]. Note also that the con- 
formal coupling ^ = — 1/6 corresponds to a model where 
there is actually no scalar degree of freedom despite the 
fact that there seems to exist one in the parametriza- 
tion ()1.2p (see [s^). The free parameter s has a physical 
meaning since it is connected to the ratio of critical points 
xi and and to the coupling. For this form of F{^) we 
can now determine A2 from ()3.46p or p.47p and B2 from 
(|3.44p and finally arrive at a solution for a and (3. For the 
two values of ^ and eqs. (|3.46p and ()3.47p we get three 
degenerate solutions for a and (3, ie they correspond to 
the same form of the potential U (see ea. p.33p ). The 
solutions are 



ai 



f3i 



a2 



6C 

a''AB{2s + 3)($ - $0)^^^+^ 



12C(s + 1) 



12C(s-M) 
a'AB{2s -I- 3)($ - $0)^^ 



+1 



12C(s-K 1) 
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Numerical results 



^3 



/33 



a"AB(2s + 3)($o-*)' 



12C(.s + 1) 

It is easy to show that for all three cases we have 

and that from eq. p.34p we get 

= [/o($ - $o)^^ = C/o($ - $0)^"^'^ (3.51) 

where Uq is a constant determining the scale of the poten- 
tial and is not directly measurable, but can be rewritten 
in terms of observable parameters like Hq, qq, fim. A 
noteworthy feature of eq. (|3.5ip is that it exactly coin- 
cides with (I3.8p , thus hinting towards a non-trivial phys- 
ical content in this class of scalar-tensor Lagrangians. 

In order to find the solution to the field equations we 
need to find the value of the constant of motion Eq from 
eq. ([3:201 . This is, 



Sn = 



fe(2s + 3)2 d(a''+3($- $0)^^^+^) 



(s + l)2(s + 2)(s + 3) 



dt 



where b = — . Integration of (|3.52p yields 



(3.52) 



$ - $0 a"~c^+5Ji2=^+8.+6 (3.53) 



where c : 



b(2s-|-3)2 ) 

Plugging (|33T|) and ([333| into jM]) we can get a(t) 
and <^{t): 



a(t) t=+3 ; — , — -T7T^5 h ao 



1+^ 



(s + 6)(2s + 3)2 



(3.54) 



and 



/ 8cs(s-H)(s + 2)(s + 3)2£/o^^ 
I (s-h6)(2s-f 3)2 



ao 



(3.55) 



where ao is an integration constant. 

Due to the structure of the above general solution, 
the cases s = 0,-1,— 3/2,— 2 and s = —3 have to be 
considered apart. The solutions for s — and s = —3 
correspond to the minimal coupling where F = Fq and 
U — A and to the quartic potential case where ~ $2 
and [/ ~ In these situations, the solutions assume 
oscillating or exponential behavior (for a discussion see 



As final step, in order to confirm the validity of our 
analysis, we perform a numerical evolution of the dy- 
namical system to compare the form of the potential U 
and couphng F of eqs. (|3.4p and (|3.8p (found also by the 
Noether Symmetry Approach), with the corresponding 
form obtained from the numerical analysis. The steps 
involved in the comparison arc the following: 



Numerically solve the dynamical system (|2.22p . 
(|2:24)) . ((2:271) and obtain xij N) X2iN), 

x'^{N) and X4{N). Integrate equation (|2.16p to get 



F{N) ^ Foer^^^.^'"''-^"^'^'^' 
and use eq. (|2.18p in the form 

$'(Ar)2 ^ Qxl{N)F{N) 



(3.56) 



(3.57) 



to obtain ^(N). The resulting form of F(<i>) in 
both the numerical (red continuous line) Eq. (|3.56p 
and its analytical approximation (blue dotted line) 
Eq. (|3.4p is shown in Fig. 4. 



• Use equation (|2.17p to obtain U{N) numerically 
(red continuous line in Fig. 5) 



U{N) = 3x2iN)F{N)H{N)'^ 



(3.58) 



and compare with the analytical form (blue dotted 
line in Fig. 5) of Eq. 



In the methodology mentioned above we used the ex- 
act coefficients xi, X2, xi and xa from Table I in each 



era for the analytic forms (|3.4p and (j3.8p . Also, the 
initial conditions used in the numerical evolution were 
F{N = -30) = ^{N = -30) = 1. As another test, 
we fitted the numerically obtained F{N) and U{N) from 
eas. (j3.56p and (|3.58p respectively to obtain the coeffi- 
cients of the analytic forms p.4p and (|3.8p and found 
that the results were in good agreement (see Figs. 4 and 
5 the green long-dashed lines). A noteworthy feature of 
Figs. 4 and 5 is a small plateau that appears during 
the deSitter era (see the zoomed region in Fig. 4). The 
reason for the existence of the plateau is that as the sys- 
tem evolves towards the deSitter era the coupling F(<i>) 
"freezes", since xi{Nds) (see eq. (|2.16p and Fig. 
2a), much faster than the field 



IV. CONCLUSIONS AND OUTLOOK 

We have reconstructed the form of the gravitational 
coupling -F($) and the potential [/(<&) of scalar-tensor 
quintessence by demanding that it reproduces a ACDM 
cosmic history through the radiation {weff — matter 
{Wf,ff = 0), and deSitter (we// = —1) eras. We have 
found that apart from the usual general relativistic solu- 
tion with a constant coupling F(<I>) = Fq and potential 
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10 20 30 

Log(O) 



FIG. 4: The form of log{F{^)) in the numerical reconstruction (red continuous hne), its analytical approximation (blue dotted 
line) and a fit of the numerical reconstruction using Ea. (|3.56|) (green long-dashed line). The agreement between the three 
approaches is very good. The reason for the existence of the small plateau, see the zoomed region, is that as the system evolves 
towards the deSitter era the potential -F("l>) "freezes" much faster than the field 




10 20 30 

Log(O) 



FIG. 5: The form of log{V{^)) in the numerical reconstruction (red continuous line), its analytical approximation (blue dotted 
line) and a fit of the numerical reconstruction (green long-dashed line). The agreement between the three approaches is very 
good. The potential exhibits a small plateau in the deSitter era for the same reason as -F($) (see caption of Fig. 4) 
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U{^) = Uq (corresponding to Newton and cosmologi- 
cal constants) , there is another consistent solution which 
reproduces the same cosmic history. According to this 
solution 

= e($-C)2 (4.1) 
[/($) ^ F($)" (4.2) 

where m, C are arbitrary constants (m however is nega- 
tive). In this new solution the 'radiation', 'matter' and 
'deSitter' expansion rates, for the 'attractor' trajectory 
shown in Figs 2 and 3, is dominated by dark gravity 
through all epochs. This is indeed a potential problem 
for this type of trajectories but it could also be a potential 
blessing since this type of solutions have the correct ex- 
pansion rate at all epochs without the use of dark matter 
or dark energy. A proper test of these models for detailed 
comparison with observations would require analysis of 
large scale structure formation (analysis of evolution of 
perturbations). Such an analysis is beyond the aims of 
our present analysis but it is an interesting extension of 
this project. 

We should also stress that not only the attractors but 
all the fixed points we found are potentially interesting 
because as discussed above a more physical stability anal- 
ysis would fix F{N) and allow H{N) to vary thus intro- 
ducing and eliminating instability modes. Thus the phys- 
ically interesting part of our analysis is the actual values 
of the fixed points and not their stability which could 
change if H{N) were allowed to vary. On the other hand, 
going back to the considerations developed in Sect. II re- 
garding the number of parameters beside m, we have to 
stress again that if we had used the parameter n ~ F $ /F 
in our analysis (thus fixing F{^) and allowing H{N) to 
vary), the autonomous system would be very different 
and, depending on the values of n, the stability and the 
form of H{N) would correspondingly vary. An alterna- 
tive approach could be to fix F and U, as in Eas. (|4.ip 
and (1321), and allow H{N) to vary about a the ACDM 
background, along the lines of Ref. [26|. Very likely, this 



approach could provide a comprehensive stability analy- 
sis of the ACDM model since it involves also the proper 
fluctuation modes of H{N). This extension is out of the 
lines of this work and will be faced in a forthcoming pa- 
per. 

Another point is that phantom behavior can be easily 
realized in scalar tensor theories, see Refs [ill, UHl for a 
discussion, and this is an attractive feature of these theo- 
ries. In fact we could have chosen to reconstruct different 
forms of H{N) giving late time phantom behavior thus 
deriving different forms of F{N) at late times with differ- 
ent fixed points. This late time reconstruction has been 
undertaken in Refs. [l3| and [s^. However, clearly the 
early times behavior of the reconstructed F{N) would 
be unchanged even in the phantom case. Fixing F{^) 
could also lead to phantom behavior but we would have 
to guess a proper form of F. 

A completely independent way that can also lead to 
the form of the new solution presented here is obtained 
by imposing maximal Noether symmetry on the scalar- 
tensor Lagrangian. We have demonstrated that impos- 
ing such a symmetry leads uniquely to exactly the same 
form of potentials as (14. 1|) and (j4.2|) . It also leads to a 
conserved charge So which allows the derivation of exact 
solutions for the evolution of the scale factor a{t) and the 
scalar field ^{t). 

This intriguing coincidence of the two approaches hints 
towards a non-trivial physical content in this class of 
scalar-tensor Lagrangians. It is therefore important to 
study the evolution of cosmological perturbations in this 
class of models in order to test them using Large Scale 
Structure and CMB observations. 
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